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Abstract: In this paper we present a mathematical analysis for a steady-state laminar bound- 
ary layer flow, governed by the Ostwald-de Wael power-law model of an incompressible non- 
Newtonian fluid past a semi-infinite power-law stretched flat plate with uniform free stream 
velocity. A generalization of the usual Blasius similarity transformation is used to find similar- 
ity solutions m. Under appropriate assumptions, partial differential equations are transformed 
into an autonomous third-order nonlinear degenerate ordinary differential equation with bound- 
ary conditions. Using a shooting method, we establish the existence of an infinite number of 
global unbounded solutions. The asymptotic behavior is also discussed. Some properties of 
those solutions depend on the viscosity power-law index. 
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1 Introduction 

In view of their wide applications in different industrial processes, and also by the interesting mathematical 
features presented their equations, boundary-layer flows of non-Newtonian fluids have motivated researchers 
in many branches of engineering in recent years. The most frequently used model in non-Newtonian fluid 
mechanics is the Ostwald-de Wael model (with a power-law rheology ll2l[3l|4l|5l|6l), which the relationship 
between the shear stress and the strain rate is given as follows 



'xy 



k\Uy\'' ^Uy, (1) 



for n = 1 the fluid is called Newtonian with dynamic coefficient of viscosity k. For n > 1 the behavior of 
the fluid is dilatant or shear-thickening and for < n < 1 the behavior is shear- thinning, in these cases the 
fluid is non-Newtonian and k is the fluid consistency. In this work we shall restrict our study to the dilatant 
fluids, then throughout all the paper, the exponent n will be taken in the range (1, oo). The problem of lam- 
inar flows of power-law non-Newtonian fluids have been studied by several authors. For the sake of brevity, 
we mention here some examples, Acrivos et al.Q and Pakdemirli lEl derived the boundary layer equations 
of power-fluids, Mansutti and Rajagobal |l9| investigated the boundary layer flow of dilatant fluids. Adopt- 
ing the Crocco variable formulation, Nachman and Talliafero 1101 established existence and uniqueness of 
similarity solution for a mass transfer problem. Filipussi et al. ifTTI obtained similarity solutions and their 
properties using a phase-plane formalism. Recently numerical solutions have been given by Fee and Biiyuk 
in IIT2I for the steady laminar free convection over a heated flat plate. 

More recently Guedda iTTS] studied the free convection problem of a Newtonian fluid, he showed the exis- 
tence of an infinite number of solution and studied their asymptotic behavior. In this work we aim to extend 
the analysis of |[T3l to the non-Newtonian case, we are interested also in the effect of the power-law index 
on the existence and the asymptotic behavior of solutions. 

The remainder of this work is organized as follows, in the next section, we introduce the mathematical for- 
mulation of the problem, section 3 deals with some preliminary tools which will be useful in section 4 and 
5 to prove the main results. Finally, we give some concluding remarks in section 6. 



2 Similarity procedure 

The problem is geometrically defined by a semi-infinite power-law stretched rigid plate, over which flows a 
non-Newtonian fluid obeying to ©. The main hypotheses for the mathematical formulation of this problem 
are given by: 
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• Two-dimensional, incompressible and steady-state laminar flow, 

• Physical properties ai^e taken as constants, 

• Body force, external gradients pressure and viscous dissipation are neglected. 

Under these assumptions, and referred to a Cartesian system of coordinates Oxy, where y = is the plate, 
the X-axis is directed upwards to the plate and the y-axis is normal to it, the continuity and momentum equa- 
tions can be simplified, within the range of validity of the Boussinesq approximation Q, to the following 
equations 



UUx + VUy = v{\UyY' ^Uy)y, 



< 



(2) 



Ux+Vy = 0, 



The functions u and v are the velocity components in the x— and y— directions respectively. 
The boundary conditions accompanied equation Q are given by 

u{x,0) = Uw{x), v{x,0) = Vu,{x), u{x,y)^{) as y — > oo. (3) 

The functions : C/^(x) = u^x"^ is called the stretching velocity and Uw > 0, the exponent m is negative, 

m{2n— 1) — n 

and Vw{x) = v^x is the suction/injection velocity where Vw > for suction and < for 

injection. 

From the incompressibility of the fluid we introduce the dimensionless stream function ^/J = y) satis- 
fying (n = 'lpy, V = -Ipx)- 

Hence equations Q are reduced to the single equation 

'4>y'>Pxy - -Ipxtpyy = '^{\ijyy\''~^'^yy)y (4) 

The boundary conditions © are transformed into 

m(2n — 1) — n 

■0y(x,O) = ipxix,0) = -v^x ^^+1 , ?/;j^(x,y) ^ as y ^ oo. (5) 

Since the broad goal of this paper is to obtain similarity solutions to (H)),© we introduce the following 
similarity transformations 

ij{x,y) ■.= Ax"f{t), t:=By^. (6) 

xi^ 

Where A,B,a and (3 are real numbers, / is the transformed dimensionless stream function and t is the 
similarity variable. In terms of equation (|4| can satisfy the ordinary differential equation of the shape 
function : / 

(irrVT +«//" = («-/?)/'', (7) 

where the primes denote differentiation with respect to t, if and only if the following 

a{2 — n) + f3{2n — 1) = 1, and a — /3 = m, 
holds, and the parameters A, B and v satisfy 

^^n-2^2n-i ^ and m = a-(3. (8) 

It follows that 

1 + m(2n — 1) ^ 1 + m(n — 2) m(2n — 1) — n 

a = — i , P = — i , and p= — . 

n+1 n+1 n+1 

Consequently, we have 

1 l + m(2n — 1) 1 l + T?i(n-2) 

'ip{x,y) := 1^"+^ X /(^)) t:=u~"+^yx . (9) 
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The corresponding boundary conditions ^ are expressed as 

/(0) = ^, /'(0) = ^, /'(oo) = lim m = 0. (10) 
In the remainder, we deal with the following problem 

(irrv")'+«//"- w = o, 

(11) 

/(O) = a, /' (0) = b, /'(oo) = limt^oo f'{t) = 0. 
For Newtonian fluids (?i = 1), problem (ITTI ) reads 

f f"' + aff"-mf" = 0, 

(12) 

t /(O) = a, /'(O) = 6, /'(oo) = 0. 

We notice that this problem arises from two different contexts in fluid mechanics when looking for similarity 
solutions. First, in natural convection along a vertical heated flat plate, embedded in a saturated porous 
medium, where the temperature is a power function with the exponent m, for more details, we refer the 
reader to |[T3l [141 [161 and the references therein. Equation ([T2l ) appears also in the study of the boundary 
layer flow, of a Newtonian fluid, adjacent to a stretching surface with a power-law velocity (see |[T5l [TTl ). 
In |[T4l . the authors proved that ([T2l) with a = b = 0, has a solution (which is bounded) if m > — | and 
this solution is unique for < m < g. In |[T5l the author gives a complete study about existence and 
nonexistence of solutions to ([T2l) . where 6 = 1. 

Recently some new results have been obtained in |[T3l . The author considered the problem ([T2l) . where 
m G {—a, 0). He showed that, under some assumptions, problem ([T2l) has an infinite number of unbounded 

a 

solutions and these solutions satisfy f{t) ~ t^-m ^ as t goes to infinity. 

Based on the results of [il3il . the interest in this work will be in existence and asymptotic behavior of solutions 
of problem ([TTl ). 

3 Preliminary results 

As it is announced above, the existence of solutions will be established by a shooting method. We replace 
the boundary condition at infinity by /"(O) = d, where d 0. Therefore, we consider the initial value 
problem 

(irrv")' + «//"-"^/'' = o, 

(13) 

f{0) = a, f'{0)=b, f"{0) = d. 

We shall see that for appropriate d problem ([T3] ) has a global unbounded solution and this solution satisfies 
the boundary condition at infinity. 

Remark 1 We notice that for n ^ 1, equation ([TTI) i can be degenerate or singular at the point to where 
/"(to) = 0. The existence the to is done for d > 0. We shall see also that /'" is not bounded at to (the 
solution / is then not classical). By a solution to ([TTTl we will mean a function / G C^(0, cxo) such that 
|/"|"~^/" G C^(0, oo), /'(oo) = and /"(oo) = 0. Note also that any solution is classical on any interval 
where the second derivative does not change the sign. 

Consider now the initial value problem ([T3] ) with n > 1, a, d G M, 6 > and ?n G {—a, 0). 
By the classical theory of ordinary differential equations the above problem has local (maximal) solutions 
on some interval (0, T^), < oo and they ai^e uniquely determined hy d {d ^ 0). Let us denote this such 
solution by fd- Integrating ([T3l) i to het the following identity 

l/d r"Vd (t) + afd{t)fd{t) = \dr-'d + aab + (m + a) f f'd{sfds, ^ t < Td- (14) 

Jo 



which will be used later for proving some results. 

A solution fd of ([IB, is of class on [0, T^), and satisfies |/^T'~Vd e C^{[0, Td))- We shall investigate 
whether fd admits an entire extension. First, we give the following result characterizing the existence time 
Td. 

Proposition 2 Let fd be the local solution to M3^ , if Td is finite then the functions fd, f'd and f'd are 
unbounded as t approaches Tdfrom below. 

Proof. Similar to lIBlfTSl. ■ 

Let us note also that if we require a classical solution of (fT3l) (ie. / G C'^(0, oo)), it is possible that / ceases 
to exist at some T < oo and such that /, /' and /" remain bounded on [0, T). More precisely we have the 
following result. 

Proposition 3 Let fd be the local solution to H13i where n > 1 ans d ^ 0. Assume that there exists 
to G (0, Td) such that fjito) = 0. Then d > 0, < on (to, Td) and f'J^' is unbounded on (0, to)- 

Proof. Assume first that d < 0. Therefore < on [0, e), e small, and the following equation 

ni/^rVd+«/r- w = o, (15) 

holds on (0, e). Hence 

{fyy ='^e^\f'l\'-V',\ on (0,e), (16) 

where 

F{t) = - f fd\f'l\'~''{s)ds. 
n Jo 

Consequently, the function t f'J^e^{t) is decreasing, and then /^'(t) remains negative for all t € [0, Td). 
A contradiction. Then d > 0. Actually, we have > 0, > 6 on (0, to) and equation (ITSl) holds on 
(0, to). Assume now that is bounded on (0, to). Thanks to equation dUJi we deduce that /^(to) = this 
is contradiction with /'(O) > 6. ■ 

From the above we can see, in particular, that < on (0, Td) for any d < 0. Then fd G C°°([0, Td)). 
While for the case d > the solution fd is not classical. 

Proposition 4 Let fd be the local solution to i!d3}for d and n > 1. If Td < oo then limt^j-^ fdit) = 
—oo. 

Proof. First we show that supp^y^) \fd{t)\ = oo. Suppose not and /^'(to) = holds, for some to G (0, Td). 
From (fT4l ) we get 

-{-fd)''{t) + afd{t)fd{t) = af'd{to)fd{to) + {m + a) f f'd{sfds, y to < t < Td. 

J to 

Hence 

^/J(t)-a/^(to)/,(to)(t-to)-^/|(to) = (m + a) /* f f'd\s)dsdT + [\-f'J)^{s)ds. 

^ ^ Jto Jto Jto 

Since the right-hand side of the above is monotonic increasing with respect to t, the function fd has a finite 
limit as t ^ Td. Consequently the function {—fd )"' is integrable on (to, Td). Since n > 1 we deduce that 
f'J^ is also integrable on (to, Td). Therefore the function is bounded. Next we use ([T4l ) to deduce that /J 
is also bounded. A contradiction with Proposition |2l 

It remains to prove that the hypothesis f'^ > on (0, Tf^) leads also to a contradiction. Actually, in such 
situation, we know that fd is classical and satisfies (fTSl) . which yields to 

/ fll\n-2 fill ^ ^ f 
[Id) Id < ~-Jd, 
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and then 

(/^r-V7<- sup 

Therefore and are bounded. A contradiction. 

Because /rf is monotonic on (r, T^), for some < r < T^, we deduce that goes to infinity as t — > T^. 

Finally, to show that fd{t) — oo a.st Td we assume on the contrary that limf^j-^ fd{t) = oo. Hence the 
functions fd and are positive on (r, T^). Moreover, using (fTTTl i . we can deduce from the Energy-function 
defined by 

E{t) = ^\mr'-^f:i', (17) 

n + 1 3 

and satisfies E'{t) = —afdfd'^ < 0. That and ai^e bounded. Hence /^^ is also bounded and this is a 
contradiction with Proposition^ We conclude that if Td is finite the function fd goes to minus infinity as 
t — > oo. ■ 



4 Existence of solutions 

In this section we shall obtain a sufficient condition on d such that the local solution fd of (fT3] ) is global and 
satisfies the condition /^(oo) = 0. We show that, for each d satisfying \d\'^~^d > —aab, fd exists on the 
entire positive axis and satisfies /^(oo) = 0. We begin by a simple observation that: if m + a > and 
\d\'^~^d > —aab, (fT4l) yields the important fact that fd cannot have a local maximum. Thus we prove the 
following result. 

Theorem 5 Let a € M, 6 > and m G {—a, 0). For any d such that jdj'^^^d > —aab, there exists a unique 
global solution fd, to U3i . which goes to infinity with t, and its first and second derivative tend to Q as t 
approaches infinity. 

For our analysis, we need to distinguish two cases for the parameter a = fd{0)', namely a > and a < 0. 
First we prove the following lemma. 

Lemma 6 If a > and \d\^~^d > —aab the functions f'^ and fd are positive on (0, Td) and Td = oo; that 
is fd is global. Moreover f'^ and f'J are bounded. 

Proof. Because \d\"'~^d-\-aab > 0, the first assertion of the lemma is immediate from (fT4l) . To demonstrate 
that Td = oo it suffices to show that fd remains bounded on any bounded interval [0, T]. Let us consider the 
Lyapunov function E for fd defined by ([TT] ). Since 

E'it) = -afdff < 0, 

thanks to (fTTTl i . it is seen that 

^^\mr' - jm' < ^,\dr' - w < t,. 

This in turn implies that f^, f^ and then fd are bounded on [0, T]. ■ 

Lemma 7 If a > and \d\"'~^d > —aab, fd{t) tends to infinity with t, f^ and f'^ tend to zero as t ^ oo. 
Proof. Since is monotonic on {ti, oo), ti, large enough, and bounded there exists a Z > such that 

lim /^(t) = I. 

t—tOO 

This implies the existence of a sequence tending to infinity with n satisfying lim„_>+oo f'di^n) = and 
then lim^^oo (0 = 0) with the help of the energy function E. 
Now we assume that fd is bounded, therefore / = 0. Subsequently 

;>oo 

\d\'^-^d + aab+{m + a) j f'^{tfdt = d. 

Jo 



This is impossible. Therefore fd is unbounded and then limj_^+oo fd{t) = oo. It remains to prove that I = 0. 
Assume on the contrary that ^ > 0. Together with (fT4l ) we get 

l/dT"Vd W = + + a)lH + oit), 

\f'ir'f'M=mlh + o{t), 

as t approaches infinity, that This is only possible if m = 0. Consequently I = 0. m 
Next we consider the case a < 0. 

The first simple consequence is that fd{t) < and /^(t) > for small t > 0. Since fd cannot have a local 
maximum, we have two possibilities 

• Either fd{t) vanishes at a some point and remains positive after this point. 

• Or fd{t) < Oforallt > 0. 

Hence the proof of Theorem [5] is completed by the following lemma. 

Lemma 8 Assume a < and \d\'^^^d > —aab. Then fd has exactly one zero, goes to oo with t, and the 
functions f'^, f^ converge to as t ^ oo. 

Proof. Assume that the first assertion holds. Since is positive we deduce that fd is bounded and then is 
global. On the other hand, using (fT4l ) one sees that f'J > 0. Therefore we get lim(_>oo fd{t) ^ (ojO] and 
limi^oo fdi^) = 0, since is monotonic. This is absurd since is positive and increasing function. Hence 
fd has exactly one zero, say to- To finish the Proof of Lemma [8] and therewith that of Theorem [5] we note 
that the new function 

Ht) = fdit + to) 

satisfies equation (fTT]) i and 

h{0) > 0, h"{0) > -ah{0)h'{0). 
Therefore we use Lemmas [6] and |7] to conclude. ■ 

In the next result we complete our analysis on the existence of global solutions by the case 5 < 0. 
Theorem 9 Let 6 < 0, a > and m G (—a, 0). For any d > satisfying 

ad"" - hf^d""-^ + aaH > 0. (18) 

The unique local solution, fd to di is global unbounded and satisfies lim^^oo f'di'^) = linij^oo f'di't) = 0. 

Proof. Since a, d > and 6 < 0, there exists a real to > such that fd is positive, decreasing and convex 
on (0, to) - Define 

T = sup [t : fd{s) > 0, f'dis) < 0, f'Jis) > 0, for all s G (0, t)} . 
The real number T is larger than to and may be infinite. 

Assume that T = oo. Then the function fd has a finite limit at infinity and /^(t) (and /^') go to zero as 
t —>■ oo. Since the function 

H = fdif'ir'f'j - \f'd\f'ir' + o^fdfd, 



satisfies 



H' = fif 



l\2 



a(n — 1) 

m + 2a-\ ^ 

2n 



m{n — 1) 
2n 



ifdnfir' 



thanks to ([TTIli . we deduce that H is increasing on (0, oo). Hence for t > we have 

H{0) < lim H{t) = 0, 

i— »oo 

which yields to 

ad^ - ife^rf""^ + aa^b < 0. 
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A contradiction. Therefore T is infinite. Next, we assume that fd{T) = or f^{T) = 0. Arguing as above 
we deduce H{0) < and then we get a contradiction. In conclusion if condition ( fTSl ) holds the function 
fd has a local positive minimum at some ti > 0. We use Theorem [5] to deduce that the new function 
h{t) = fd{t + ti) is global, unbounded and satisfies h'{oo) = h"{oo) = 0. The proof is finished. ■ 



Remark 10 We notice that we can extend the results of Theorem [5] to the case (—2a, —a) by using the 
function H defined above, as in the work by Guedda |[T3ll for the Newtonian case. 



5 Asymptotic behavior 

In this section we shall derive the asymptotic behavior of any possible global unbounded solution to (fTTI ) for 
m € (—2a, 0). First we give the following result 

Lemma 11 Let f be a positive solution to Alii for m G (—2a, 0). Then f goes to zero at infinity and f" is 
negative. 

Proof. Since / is monotonic on [to, oo), Iq large enough, we get the positivity of /' and / on (to, oo). 
In addition we use the Lyapunov function to get the boundedness of /' and /". Arguing as in the previous 
section we get that /' — > and /" < for large t. m 

Proposition 12 Assume that n > 1 and m G (—2a, 0). Let f be a positive solution to diil) . Then 

lim /,(t)/^(t) = lim (|/"r-i/")'(i)=0. 

I— ►oo t— >oo ' 

Proof. Thanks to lemma ([TB we have f'{t) > 0, f"{t) < for all t > to, to large enough and /' and /" 
tend to as t — 5- cxD. Then equation ([TTI li can be written as 



r + -/ri/ 



n 



n 



By differentiation we have 

j{iv) _|_ Jill 



a(2 - n) |j//|l-nj _ "^(1 -"-) |j//|-n-lj//j/2' 



n 



n 



a -2m „,„,,„,x-n 



n 



-rnr 



(19) 



Then the function f"'e is monotonic increasing on {to, oo), where 

^, ^ a(2--n) _ ?7i(l - n) „i^„^,2 

n n 

This indicates that the function /'" has at most one zero. Because /" is negative and goes to at infinity, we 
deduce that f"'{t) > on (ti, oo), for ti large. On the other hand, from (fTTb i we deduce 



(|/'r-VT + («-2m)/7" 



-a//" 



(l/"r ^f"y-\- " f'^ is positive and monotonic decreasing on (inf {to, ti} , 00) 



Therefore the function t 
Together with the fact that /' tends to as t ^ 00 we deduce that 

lim (irr-v")'(t) = o 

t— »+oo 

and then we conclude that ff"{t) ^ as t — > 00, thanks to (fTTb i . ■ 
Proposition 13 Let f be a solution to f l77] ) where n > 1, m G (— 2a, 0). Then 

00, if m + a > 0, 



lim fr{t)={ 



L G (0,00) 
I 0, 



if m + a = 0, 
if m + a < 0. 



Proof. Let / be a global solution to ([TT]) . First we claim that there exists tQ > such that 

\f"r'nto)+aff{to)>o. 

Suppose not; that is 

irrV'w +«//'(*) <o, 

for all t > 0. Since f"{t) ^ as t ^ oo the following 

/" + aff < 

holds on some (ti, oo), ti large. Consequently the function /' + ^/^ is decreasing and goes to infinity with 
t, which is absurd. Now we use the identity 

+ affit) = irr'nto) + aff (to) + (m + «) /* f'\s)ds, 

to deduce that //' has a limit L € [0, oo] at infinity. This limit is finite for a + m = 0. Assume that 

a + m 7^ 0. If L G (0, oo) we get immediately that /' ~ at infinity which implies that //' —i- oo. 

A contradiction. Consequently L € {0, oo}, we use again the above identity to conclude that L = oo if 
m + a > and L = Oifm + a<0. ■ 



Remark 14 We stress that the condition 

irrV"(io)+«//'(to)>o, /'(to)>o 

is necessary and sufficient to obtain a global solution converging to plus infinity with t in the case m G 
i-a,0). 

Now we are ready to give the result concerning the large t— behavior of solutions to (fTTI) . 

Theorem 15 Suppose n > 1, —2a < m < 0. Let f be a solution to diil ) such that / — > oo. Then there 
exists a constant, A > 0, such that 

f{t)=t^{A + o{l)), (20) 

as t ^ oo. 

Proof. Let / be a global solution to ([TT]) . First we prove the result for the case m + a > 0. 

Let to be a real number such that /" < and > on (to, oo). Dividing equation (fTTTl i by //' gives 

{\f"r~'f"y f f" 

= m a — . 

//' / /' 

Integrating over (ti, t), for ti > to, gives 

'{\f"r-'f"y 



ff 



-ds = log (r (t)/'""(t) - log (/™(ti) - /'""(ti)). 



According to Proposition [131 //' goes to infinity with t and then the left hand side of the above is integrable, 
consequently has a positive finite limit at infinity. The desired asymptotic behavior ( [201 ) follows by 

a simple integration. 

Now we deal with the case m + a < 0. For this sake we define 

* = ^(/)i/'r-V" - \v'{m'?\f"r^ + a^iDif, 

where is a smooth function. Then if follows from ([TTl li 

V{t) = f [af^'U) + (a + m)^] - ic^"(/)/' V'T"' " ^^'(/)/' V'T-V"/'" ■ 
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Let the function if be defined by 



ip{S) = S a . 



It satisfies the following differential equation 

asip'{s) + (a + m)ip = 0. 

This implies that 



*' = -2^"(/)/ 1/ 



2af 



■ifT\f"r' + aff' 



and then 



f/\n-2 fill 



for all t > to, to large. Therefore is integrable on [0, oo) and then has a finite limit at infinity, say L. 
Next we show that L > 0. It will be sufficient to show that ^'(ii) > for some ti large. Suppose not; that 
is for any t > t2, t2 large we have 

irr-V" - ^^/'V'T"' + «//' < 0. 

2a/ 



Since a + m < then 



/" + aff < 0, 



from which we deduce, as above, that /' + ^/^ is a decreasing function going to infinity with t. A contra- 



2 • 

diction. We conclude that limt^oo f~~ f = Finally, a simple integration leads to estimate (l20l) . 
To finish, we pay attention to the case m = —a. Here the identity (fT4l) leads to 

i/"rV"+a//' = i7rs+«a&. 

According to Theorem |5l / is global and satisfies / ~ ta at infinity. ■ 



6 Conclusion 

The laminar two-dimensional steady boundary layer flow, of a non-Newtonian incompressible fluid, over 
a stretching surface have been considered. Using the shooting method, existence of global unbounded 
similarity solutions have been shown, the dependency of those solutions on the power-law index has been 
investigated, and their asymptotic behavior was also discussed. 

Coming back to the original problem (II]),© we find that, for —2a < m < 0, the stream function satisfies 

l + m(2n-l) (2-n)m-l 

^(x, y) ~ y as yx "+i ^ oo. 
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